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ABSTRACT 


The  traction  developed  by  a  wheel  can  be  determined  by  ap¬ 
propriate  integration  of  the  soil-wheel  interface  stresses.  In 
soft  soils  these  are  governed  by  failure  conditions  in  the  soil 
brought  about  by  the  wheel  load.  The  geometry  of  the  two  failure 
zones  (front  and  rear),  as  well  as  the  associated  stresses,  con¬ 
stitute  the  solution  of  the  differential  equations  of  plasticity 
for  the  boundary  conditions  defined  by  the  geometry  of  soil-wheel 
interface  and  the  interface  friction. 

The  numerical  solution  of  the  differential  equations  of  plas¬ 
ticity,  as  well  as  the  selection  of  entry  and  exit  angles  that 
yield  matching  interface  stresses  at  the  common  point  of  the  front 
and  rear  failure  zones,  are  amenable  to  the  application  of  com¬ 
puter  techniques.  A  flow  diagram  shows  the  computation  of  fail¬ 
ure  zones  and  associated  interface  stresses  that  yield  the  load, 
torque,  and  traction  for  a  given  wheel  geometry  and  interface 
■friction.  Slip-shear  stress  equations  are  used  to  relate  inter¬ 
face  friction  to  slip. 

The  solution  of  the  practical  problem  of  determining  traction 
and  slip  for  given  soil  conditions  and  applied  wheel  load  and  tor¬ 
que  requires  the  iteration  procedure.  This  is  shown  in  another  flow 
diagram.  Inputs  for  the  computer  program  schematically  shown  in 
the  flow  diagram  are  soil  properties,  wheel  load,  and  torque;  the 
end  outputs  are  traction  and  slip,  while  entry  and  exit  angles, 
interface  stresses,  and  geometry  of  the  failure  zones  can  be  ob¬ 
tained,  if  desired,  either  in  the  form  of  printed  output  or  data 
files  suitable  for  viewing  on  a  visual  display  terminal. 

For  the  computation  of  traction  exerted  by  tires,  a  soil- 
tire  model  has  been  developed  that  allows  for  the  deflection  of 
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the  tire  and  attendant  restraint  of  the  normal  stresses  at  the  in¬ 
terface.  The  front  and  rear  failure  zones  are  separated  by  a  zone 
in  which  the  soil  is  not  in  failure  condition  and  the  normal 
stress  is  the  maximum  compatible  with  the  inflation  pressure  and 
carcass  strength. 

Traction  values  computed  on  the  basis  of  soil  failure  condi¬ 
tions  by  the  above  computer  program  are  compared  with  experimental 
results.  The  agreement  is  reasonably  good  when  soil  strength 
properties  determined  by  triaxial  tests  are  used. 
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“  cohesion 

«=  drawbar  pull  or  drag 
=  slip 

*  constant  defining  threshold  slip 
=  constant  in  slip  equation 

*  load 

«  limit  normal  stress  on  tire-soil  interface 
=  normal  stress  measured  in  forward  field  at  a 

tn 

=  normal  stress  measured  in  rear  field  at  a 

m 

=  torque 

**  geometric  coordinates 
=  central  angle 

=  angles  defining  start  and  end  of  tire  deflection 

=  angles  defining  ends  of  flat  portion  of  tire 
deflection 

«*  entry  angle 
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«  angle  of  separation 
■  rear  angle 

y  =  unit  weight  of  soil 

6  =  angle  of  inclination  of  resultant  stress  to  normal 

e  =  slope  angle 

A  =  computed  load/input  load  (zero  subscript  indicates 

initial  value) 

<p  *  friction  angle 

■tp  =  c  cot  4> 

6  *»  angle  enclosed  by  major  principal  stress  and  x  axis 

lx  b  45°  -  <p/2 

OpC^ >°3  e  principal  stresses 

a  -  (a^  +  f 

cd  <=  DB/L  (zero  subscript  indicates  initial  value) 

1.  INTRODUCTION 

The  prediction  and  analysis  of  the  tractive  performance  of  rigid 
wheels  on  soil  is  one  of  the  basic  problems  in  the  field  of  soil- 
vehicle  mechanics.  The  soil-wheel  interaction  theory  presented  in 
this  paper  relates  traction  performance  directly  to  the  strength 
properties  of  soil,  providing  thereby  the  r-.1vtical  tool  essential 
for  a  predictive  method.  Presentation  of  the  theory  and  its  applica¬ 
tions  are  arranged  in  this  paper  in  the  following  order.  In  Section  2 
the  theoretical  background  is  presented  together  with  the  application 
of  the  theory  to  wheel  performance  calculations.  Section  3  is  allo¬ 
cated  to  the  discussion  of  numerical  solution  methods.  In  Section  4, 
results  of  experiments  are  compared  with  theoretical  predictions.  In 
Section  5,  results  of  research  work  on  rigid  wheel-soil  interaction 


are  summarized  and  further  development  of  the  theory  for  application 
to  tire-soil  interaction  is  outlined. 

2.  THEORETICAL  BACKGROUND 

The  theory  of  soil-wheel  interaction  developed  at  Grumman  has 
been  described  in  detail  elsewhere  (Ref.  1).  The  basic  concepts 
of  this  theory  are  summarized  below  for  the  convenience  of  refer¬ 
ence. 

•  Calculation  of  the  tractive  performance  of  wheels  can 
always  be  reduced  to  the  determination  of  interface 
stresses.  Once  these  are  known,  the  wheel  can  be  con¬ 
sidered  as  a  free  body  and  performance  parameters  can 
be  determined  as  resultants  of  stresses  on  a  free  body. 

•  Soil  failure  controls  the  interface  stresses  beneath  a 
driven  rigid  wheel.  Generally,  two  soil  failure  zones 
are  formed  beneath  wheels,  a  forward  one  and  a  backward 
one  (Fig.  1).  The  adjectives  indicate  both  the  loca¬ 
tion  of  these  zones  relative  to  the  wheel  and  the  di¬ 
rection  in  which  the  soil  in  the  failure  zones  tends 

to  move. 

•  The  geometry  of  soil  failure  zones  (also  called  slip 
line  fields)  and  the  associated  stresses  can  be  deter¬ 
mined  by  the  theory  of  plasticity  for  soils. 

Interaction  between  wheel  and  soil  comes  into  play  through 
the  variation  of  entry  and  exit  angles  and  through  the  effect  of 
interface  shear  stresses  on  soil  reaction.  Interface  shear 
stresses,  directly  related  to  the  applied  torque,  influence  the 
geometry  of  failure  zones  and  associated  stresses  profoundly. 

The  forward  and  backward  slip  line  fields  are  not  indepen¬ 
dent  of  each  other.  At  their  point  of  junction,  the  interface 
normal  and  shear  stresses  must  be  the  same  for  both  fields. 
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For  the  mathematical  formulation  of  the  problem,  the  follow¬ 
ing  assumptions  are  made: 

•  The  soil  is  assumed  to  be  semi-infinite,  homogeneous, 
and  isotropic. 

•  Strength  properties  of  the  soil  are  characterized  by 
the  cohesion  and  internal  angle  of  friction. 

•  Stresses  are  assumed  to  be  the  same  in  all  planes  paral¬ 
lel  to  the  plane  of  motion  of  the  wheel  (the  problem  is 
assumed  to  be  two  dimensional). 

•  Wheel  velocity  is  constant. 

•  Soil  inertia  forces  are  negligible. 

•  Soil  failure  is  assumed  to  be  governed  by  the  Mohr- 
Coulomb  yield  (failure)  criterion  as  follows: 

t  -  c  +  a  tan  (p  .  (1) 

n 

•  Pore  water  pressures  are  negligible. 

For  these  assumptions,  the  differential  equations  of  plas¬ 
ticity  for  soils,  as  derived  by  Sokolovskii  (Ref.  2),  are  valid 
and  express  the  geometric  and  stress  conditions  that  failure  zones 
must  meet.  These  equations  are 

dz  «=  dx  tan(0  ±  (i) 

(2) 

y 

da  ±  2a  tan  cpd8  *=  ~s— -  [sin(e  ±  cp)dx  +  cos(e  ±  <p)dz] 

For  properly  defined  boundary  conditions  the  above  sets  of 
the  differential  equations  of  plasticity  for  soils  yield  a  unique 
solution  in  the  form  of  slip  line  fields  and  associated  stresses. 
Slip  lines  represent  the  orientation  of  failure  surfaces  along 
which  the  Mohr-Coulomb  yield  criterion  is  satisfied,  i.e.,  the 
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Mohr  circle  touches  the  strength  envelope.  At  any  point  within 
the  slip  line  field  the  applicable  differential  equations  set 
forth  above  are  satisfied. 

There  are  some  fundamental  aspects  of  the  theory  that  need 
to  be  emphasized  here  before  proceeding  with  the  discussion  of 
the  application  of  the  theory  to  the  problem  of  wheel-soil  in¬ 
teraction.  These  are: 

•  The  Mohr-Coulomb  yield  criterion  as  applied  in  the 
above  equation  refers  to  effective  stresses,  i.e.,  for 
the  case  where  pore  water  pressures  are  negligible. 
(Apparent  cohesion  due  to  pore  water  tension  may  be 
considered  as  effective  stress.) 

•  Equations  (2)  are  valid  for  a  nonlinear  strength  en¬ 
velope  (Ref.  3)  and  may  be  expanded  to  include  pore 
water  pressures  (Ref.  4)  or  soil  inertia  forces 
(Ref.  5). 

•  The  Mohr-Coulomb  yield  criterion  implies  that  the  soil 
strength  expressed  by  Eq.  (1)  is  available  regardless 
of  the  volumetric  strain  that  is  associated  with  the 
stress  state  expressed  by  Eq.  (1).  In  soil-wheel  in¬ 
teraction  problems  in  soft  soils,  where  the  soil  is 
progressively  compressed  as  the  wheel  advances,  the 
volumetric  strain  in  the  soil  is  generally  large  enough 
to  mobilize  the  full  Mohr-Coulomb  strength. 

•  In  plasticity  theory,  solutions  obtained  by  integration 
of  Eqs .  (2)  are  termed  statically  admissible  solutions 
and  are  considered  lower  bound  solutions.  Kinematic 
admissibility  is  analyzed  by  constructing  velocity 
fields  for  the  slip  zones  on  the  assumption  that  the 
material  is  incompressible.  A  kinematically  admissible 
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solution  would  constitute  an  upper  bound.  For  soil- 
wheel  interaction  problems  in  soft  soils  the  assumption 
of  incompressibility  is  inappropriate  and  conclusions 
drawn  on  the  basis  of  such  an  assumption  are  inapplica¬ 
ble.  Experimental  evidence,  as  it  will  be  discussed  in 
Section  4,  shows  that  interface  stresses  predicted  by 
the  application  of  Eqs.  (2)  agree  reasonably  well  with 
measured  ones,  indicating  the  validity  of  the  solutions. 

Equations  (2),  with  properly  defined  boundary  conditions, 
yield  a  unique  solution  in  the  form  of  a  slip  line  field.  For  the 
passive  zone  (Fig.  1),  the  boundary  conditions  are  completely  de¬ 
fined  at  the  free  surface.  For  the  active  zone,  the  bouidary  con¬ 
ditions  at  the  interface  constitute  a  so-called  mixed  boundary 
value  problem  where  two  of  the  four  variables  (x,z,0,a)  in 
Eq.  (1)  have  to  be  specified.  In  the  solution  for  the  conventional 
bearing  capacity  problem,  z  and  8  are  specified  at  the  base  of 
the  bearing  plate.  In  the  soil-wheel  interaction  problem  it  is 
not  possible  to  specify  a  priori  the  value  of  any  two  of  the  four 
variables  at  the  interface.  Instead,  the  boundary  conditions  take 
the  form  of  a  relationship  between  x  and  z,  given  by  the  wheel 
geometry  and  a  relationship  between  9  and  the  angle  of  interface 
friction.  This  relationship  is 

9  *  |  +  i  (A  +  6)  -  a  (3) 

where 

,  ,  /sin  6\ 

a  -  arc  stn  lurr;  • 

The  interface  friction  angle,  6,  as  defined  in  Fig.  2,  es¬ 
tablishes  a  relationship  between  normal  and  shear  stresses  at  the 
interface  and  represents  the  degree  of  mobilization  of  the  soil 
shear  strength  at  the  interface.  The  orientation  of  the  principal 
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stresses  at  the  interface,  defined  by  the  angle  0,  can  be  com¬ 
puted  from  Eq.  (3)  for  any  given  6  angle. 


The  following  relationship  holds  between  the  shear  strength 
mobilized  at  the  interface  and  the  angle  5  (Fig.  2) 


tan  5 


mob 

o  +  i/ 
n  r 


(4) 


The  mobilized  shear  strength  is  frequently  related  to  slip. 
On  the  basis  of  direct  shear  tests,  Janosi  and  Hanamoto  (Ref.  11) 
proposed  the  following  relationship  between  mobilized  shear  and 
slip  for  tracked  vehicles: 


mob 


=  T 


max 


(5) 


For  compressible  soils,  which  are  of  primary  interest  in  off¬ 
road  locomotion,  this  equation  properly  describes  the  relationship 
between  shear  stress  and  slip.  When  this  relationship  is  applied 
to  the  wheel,  however,  it  is  useful  to  include  a  constant,  j^, 
in  the  slip  term  to  account  for  the  fact  that  there  exists  a 
threshold  perimeter  shear  at  which  movement  of  the  wheel  starts. 
Thus,  Eq.  (5)  is  modified  as  follows: 


mob  max 


U- 


-o+V/k, 


tan  5  =  tan  6  ( 1  -  e 

max\ 


) 

a+j0)/^ 


(6) 


The  concept  of  soil-wheel  interaction  as  outlined  before  and 
shown  graphically  in  Fig.  1  has  interesting  implications  regarding 
the  maximum  shear  strength  (t  )  that  can  be  mobilized  at  the 
interface.  According  to  this  concept,  the  soil  adjoining  the  in¬ 
terface  is  in  the  active  state  of  failure.  For  a  given  normal  and 
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shear  stress  at  the  interface  there  is  one  Mohr  circle  which 
represents  the  active,  and  another  one  xvbich  represents  the  pas¬ 
sive,  states  of  stresses.  Figure  3  shows  Mohr  circles  for  the 
active  and  passive  state  for  the  same  normal  stress  but  in¬ 
creasing  interface  shear  stress.  Stress  circles  are  shown  by 
/full  lines  for  the  active  state,  and  by  dashed  lines  for  the 
passive  state.  The  interface  shear  stress,  r,  is  shown  to  in¬ 
crease  with  the  interface  friction  angle  6  (x^  >  T2  >  Ti  anc* 

5^  >  ^  .  In  the  active  state,  the  center  of  Mohr  circles 

is  to  the  left  of  the  shear  stress  ordinate;  in  the  passive  state, 
it  is  to  the  right.  From  the  construction  of  the  Mohr  circles,  it 
is  obvious  that  the  maximum  shear  stress  that  can  be  mobilized  in 
the  active  state  is  the  one  corresponding  to  a  Mohr  circle  that 
has  its  center  at  (circle  3  in  Fig..  3).  If  the  shear  stress 

was  higher  than  this,  the  corresponding  Mohr  circle  would  repre¬ 
sent  a  passive  state.  Thus,  it  is  incorrect  to  assume  that  the 
full  soil  strength  can  be  mobilized  beneath  a  wheel  or  track. 
Passive  state  can  only  exist  beneath  a  wheel  or  track  if  the  soil 
is  pushed  toward  the  wneel  or  track,  an  obviously  meaningless 
situation  for  vehicle  mobility.  The  only  possible  stress  state  in 
the  soil  beneath  a  wheel  or  track  is  the  active  state  of  stresses 
and  it  follows  from  thxS  state  that  the  maximum  mobilized  shear 
stress  cannot  exceed  that  defined  by  Mohr  circle  3  in  Fig.  3.  In 
mathematical  terms 


T 

max 


(a  +  ip)  tan  6  =  (a  +  ip)  sin  <p 

n  max  v  n  r 


6  =  arc  tan  (sin  q>) 

max  v  ' 


(7) 

(8) 


The  assumption  of  interface  friction  angle  defines  6  at  the 
interface  and  together  with  the  geometry  defines  the  boundary  con¬ 
ditions  so  that  Eq.  (1)  yields  a  unique  solution  for  the  slip  line 
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field  geometry  and  associated  stresses.  For  the  soil-wheel  in¬ 
teraction  problem,  however,  the  boundary  conditions  for  both  the 
forward  and  rear,  slip  line  fields  have  to  be  defined  so  that  the 
solution  is  unique.  This  involves  the  determination  of  angle 
that  separates  the  forward  and  rear  fields.  The  following  con¬ 
siderations  apply  in  this  respect. 

For  dry  sand  it  was  found  by  Sela  (Ref.  6)  that  the  angle  of 
separation  equals  the  developed  friction  angle.  This  finding  is 
consistent  with  the  concept  that  the  forward  failure  zone  extends 
over  that  part  of  the  wheel  perimeter  where  the  component  of  the 
normal  and  shear  stresses,  (AD),  in  the  direction  of  motion  is 
negative  (i.e.,  resisting  the  motion),  and  the  backward  zone  ex¬ 
tends  over  That  part  of  the  wheel  perimeter  where  this  component 
is  positive.  Applying  this  concept  to  soils  with  cohesion  results 
in  the  following  relations 

AD  *  t  ,  cos  a  -  a  sin  a 
mob  n 

AD  «  (a  +  ib)  tan  6  cos  a  -  (a  )  sin  a  »  0  (9) 

n  n 

a  +  ib  , 

a  =  arc  tan  ( — - - —  tan  5  ) 

sa  \  a  ' 

n 

3.  TECHNIQUES  OF  NUMERICAL  SOLUTIONS 

The  calculation  of  wheel  performance  on  the  basis  of  the 
concepts  described  in  the  preceding  section  involves  the  follow¬ 
ing  steps: 

*  For  given  wheel  size,  soil  strength  properties,  and  in¬ 
terface  friction  parameter  (or  slip),  computation  of  the 
geometry  of  a  single  slip  line  field  and  the  associated 
stresses  by  numerical  methods. 
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M3.  ABSTRACT 


-  The  traction  developed  by  a  wheel  can  be  determined  by  appropriate  integration 
of  the  soil  wheel  interface  stresses.  In  soft  soils  these  are  governed  by 
failure  conditions  in  the  soil  brought  about  by  the  wheel  load.  The  geometry  of 
the  two  failure  zones  (front  and  rear),  as  well  as  the  associated  stresses,  are 
computed  by  solving  the.  ...differential  equations  of  plasticity  by  numerical  tech¬ 
niques.  A  flow  diagram  shows  the  computation  of  failure  zones  and  associated 
interface  stresses  that  yield  the  load,  torque  and  traction  for  a  given  wheel 
geometry  and  interface  friction.  Slip-shear  stress  equations  are  used  to  relate 
interface  friction  to  slip.  The  solution  of  the  practical  problem  of  determining 
traction  and  slip  for  given  soil  conditions  and  applied  wheel  load  and  torque 
requires  an  iteration  procedure  that  is  shown  in  another  flow  diagram.  Inputs 
for  the  computer  program  are  soil  properties,  wheel  load,  and  torque;  the  end 
outputs  are  traction  slip  and  sihkage.  For  the  computation  of  traction  exerted 
by  tires,  a  soil-tire  model  has  been  developed  that  allows  for  tine  deflection 
of  the  tire  and  attendant  restraint  of  the  normal  interface  stresses,  ;  j 
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Fig.  4  Computation  of  Slip  Line  Field  by  Finite  Differences 


*  Determination  of  a  matching  set  of  forward  and  rear  slip 
line  fields  and  the  corresponding  load,  torque,  and 
drawbar  pull  (or  drag)  . 

•  Inversion  procedure  to  find  torque  and  slip  for  given 
load  and  drawbar  pull  or  drawbar  pull  and  slip  for 
given  load  and  torque. 

3.1  For  the  computation  of  a  single  slip  line  field,  finite 
difference  forms  of  Eqs.  (2)  are  used.  These  forms  allow  the  com¬ 
putation  of  variables  at  a  point  (i,j)  from  their  values  at 
points  (i-l,j)  and  (i,j-l)  (Fig.  4).  The  procedure  is  described 
elsewhere  in  detail  (Refs.  7  and  8)  for  the  case  of  bearing  capaci¬ 
ty  computations.  For  the  problem  of  soil-wheel  interaction  it  is 
more  convenient  to  employ  a  different  sequence  of  operations.  In¬ 
stead  of  computing  the  variables  first  in  the  passive,  then  in  the 
radial,  and  finally  in  the  active  zone,  the  variables  are  computed 
along  the  first  "j"  line  in  all  three  zones  and  then  along  sub¬ 
sequent  "j"  lines  until  a  "j"  line  ends  at  a  point  at  the  in¬ 
terface  past  the  angle  of  separation.  A  "j"  line  is  then  inter¬ 
polated  between  this  and  the  previous  "j"  line  so  that  it  ends 
up  at  am  within  a  preassigned  limit  of  tolerance  (£).  This 
method  eliminates  the  time  consuming  trial  and  error  procedure  of 
finding  the  length  of  passive  zone  that  matches  the  arc  length  of 
the  active  zone.  The  length  of  the  passive  zone  in  this  procedure 
is  overestimated  so  that  the  last  j  line  would  overshoot  the  a 
angle. 

The  grid  in  this  procedure  is  larger  than  in  the  conventional 
one  (say  16  x  48  instead  of  10  x  30)  requiring  a  somewhat 
larger  computer  core,  but  the  computing  time  for  finding  a  slip 
’.xne  field  that  meets  the  boundary  conditions  at  a  is  much  less. 
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!NPUT:WHEEL  PARAMETERS  R,  B,  a  e,  a  m 
SOIL  PROPERTIES  C ,  %  tY 
TERRAIN  PARAMETER  * 

INTERACTION  PARAMETER  8 


Fig.  5  Flow  Diagram  for  the  Computation  of  a  Single  Slip  Line  Field 
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The  procedure  for  the  computation  of  a  single  slip  line  field 
is  described  in  the  flow  diagram  shown  in  Fig.  5.  The  interface 
friction  angle,  as  shown  in  the  flow  diagram,  is  assumed  to  be 
constant.  However,  the  program  can  also  accommodate  a  5  angle 
which  varies  along  the  interface. 

3.2  The  computation  of  a  matching  set  of  slip  line  fields 
starts  with  the  rear  slip  line  field  using  the  subroutine  outlined 
in  3.1,  but  allowing  for  appropriate  sign  changes  due  to  the  fact 
that  the  rear  field  is  a  mirror  image  of  the  front  field.  The 
reason  for  starting  the  computation  with  the  rear  field  is  two¬ 
fold:  first,  the  rear  angle  varies  within  a  narrow  range;  second¬ 
ly,  a  normal  stress  from  the  rear  field  (q^r) ,  at  a  ,  can  be 
generally  matched  by  a  forward  field.  This  is  not  true  vice- 
versa.  When  q^r  for  the  assumed  rear  angle  (a^)  and  interface 
friction  angle  (6)  is  determined,  the  front  field  is  found  by 
varying  the  entry  angle  (ag)  until  the  normal  stress  from  the 
front  field  (q^)  matches  q^.  within  the  allowed  tolerance. 

Since  the  interface  friction  angle  8  is  assumed  to  be  the  same 
for  both  the  forward  and  rear  fields,  the  shear  stress  at  a 

m 

from  both  fields  is  the  same  (within  the  tolerance)  when  q^  and 
q^  are  matched.  When  a  matching  set  of  rear  and  forward  fields 
is  found  by  this  procedure  the  load,  torque,  and  drawbar  pull  is 
computed  by  appropriate  numerical  integration  of  the  interface 
stresses.  The  slip  is  computed  from  Eq.  (6)  and  the  sinkage  from 
the  entry  angle.  The  flow  diagram  for  the  computation  of  a  match¬ 
ing  set  of  slip  line  fields  is  shown  in  Fig.  6. 

3.3  In  many  cases  the  problem  is  posed  differently  and  it  is 
necessary  to  solve  for  torque,  slip,  and  sinkage  for  given  wheel 
load,  drawbar  pull,  and  soil  properties.  In  the  solution  procedure 
outlined  in  3.2,  load,  drawbar  pull,  and  torque  are  functions  of  the 
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rear  angle  a  ,  and  interface  friction  angle  6,  as  expressed 
by  the  following  relationship 

L  -  f1(ar,B) 

DB  =  f2(ar,5)  (10) 

T  =  f3(ar,6)  . 

The  functions  f^,  f2,  f^  are,  of  course,  not  closed  form 
functions,  and  in  the  solution  procedure  outlined  in  3.2,  one  set 
of  L,  DB,  T  values  are  found  for  an  assumed  and  slip  cor¬ 

responding  to  an  assumed  5.  The  procedure  to  find  the  torque 
for  given  load  and  drawbar  pull  consists  of  finding  and  6 

which  yield  the  load  and  drawbar  pull;  once  the  matching  set  of 
slip  line  fields  for  these  conditions  is  found,  the  torque,  slip, 
and  sinkage  is  also  available  from  the  computations.  Even  if  the 
functions  f^  and  were  known,  the  solution  for  L  and  DB 

would  require  the  solution  of  a  system  of  two  nonlinear  equations. 
Since  there  is  no  generally  valid  theorem  for  the  solution  of  this 
problem  and  convergence  criteria  for  iterative  solutions  can  not 
be  established  if  the  derivatives  of  f^  and  f2  are  not  known, 
it  was  necessary  to  study  the  general  behavior  of  these  functions 
and  make  judicious  use  of  some  of  their  properties  to  devise  an 
efficient  convergent  iteration  scheme  for  the  solution  of  the 
problem. 

The  general  behavior  of  the  load  and  drawbar  pull  as  func¬ 
tions  of  and  5  is  illustrated  in  Figs.  7  and  8,  where  iso¬ 
metric  views  of  computed  L  and  DB  values  are  shown  for  the 
following  input  parameters: 


ir 


wheel  radius 
wheel  width 


1.15  ft 
0.36  ft 


Forward  field  friction  angle 
unit  weight 


38° 

100  lbs/cu  ft 


Rear  field  friction  angle  41° 

unit  weight  110  lbs/cu  ft 

The  problem  of  finding  the  a  and  6  values  for  given  L 
(250  lbs)  and  DB  (20  lbs)  is  also  illustrated  in  these  figures 
The  contour  line  marked  L  *  250  lbs  in  Fig.  7  represents  all 
combinations  of  the  ar  and  8  values  that  result  in  that  load, 
while  the  contour  line  marked  DB  ■»  20  lbs  in  Fig.  8  represent 
all  combinations  of  ar  and  6  that  yield  20  lbs  drawbar  pull 
The  intersection  of  this  latter  contour  line  (replotted  in  Fig.  7 
by  dash-dot  line)  with  the  L  =  259  lbs  contour  line  yields  the 
unique  solution. 


The  surfaces  representing  the  variation  of  L  and  DB  with 
ar  and  6,  shown  in  Figs.  7  and  8,  indicate  that  iteration  pro¬ 
cedures  based  on  linear  approximations  are  not  likely  to  succeed. 
Further  study  of  the  relationships  between  these  variables  in¬ 
dicated,  however,  that,  in  general,  L  monotonically  increases 
with  a  ,  and  the  ratio  DB/L  (pull  coefficient)  monotonically 
increases  with  5.  These  features  allowed  the  development  of  the 
iteration  scheme  shown  in  Fig.  9.  This  iteration  scheme  is  the 
basis  of  the  computer  program,  developed  for  the  U.S.  ATAC,  which 
computes  torque,  slip,  and  sinkage  for  given  soil  properties, 
wheel  load,  and  drawbar  pull.  The  program  requires  a  9K  core 
and  average  computer  time  is  in  the  range  of  30-60  seconds  on 
a  PDP  10  computer. 
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4.  EXPERIMENTAL  RESULTS 

In  the  development  of  the  theory,  utmost  attention  was  paid 
to  the  experimental  information  that  was  available  at  that  time 
on  the  stress  state  in  the  soil  and  at  the  soil-wheel  interface. 
Comparative  discussions  of  these  experimental  results  with  the 
theory  are  given  in  Refs.  1,  3,  and  4. 

In  the  spring  and  summer  of  1972  a  validation  test  series, 
sponsored  by  the  U.S.  ATAC  was  conducted  at  the  Stevens  Institute 
of  Technology.  In  this  test  series  a  1.15  ft  diameter,  0.36  ft 
wide  wheel  was  instrumented  to  measure  the  interface  normal  and 
shear  stresses  at  four  locations,  evenly  distributed  between  the 
centerline  and  edge  of  the  wheel.  The  tests  were  performed  in 
both  loose  and  dense  sand  and  in  a  light  clay  (loam)  at  various 
loads  and  slips.  The  strength  properties  of  the  soil  were  deter¬ 
mined  by  triaxial  tests;  cone  penetrometer  tests  were  performed 
before  and  after  the  wheel  tests  to  obtain  information  on  the 
actual  density  of  the  soil  bed.  Typical  results  of  these  experi¬ 
ments  are  given  below;  detailed  results  of  the  entire  program  will 
be  presented  in  the  final  report  on  U.S.  Army  Tank  Automotive  Com¬ 
mand,  Contract  No.  DAAE  07-72<C-0033 . 

4.1  Loose  sand.  The  friction  angle  of  the  loose  sand,  has 
been  determined  by  triaxial  tests  as  $  =  38° .  In  the  rear  field, 
however,  cone  penetrometer  tests  indicated  a  considerable  increase 
in  the  density.  The  friction  angle  corresponding  to  the  increased 
density  was  found  to  be  <t>  «  41°  . 

Fig.  10  shows  both  the  measured  and  computed  distribution  of 
average  normal  and  shear  stresses  at  the  interface  for  33  percent 
slip;  Fig.  11  for  22  percent  slip,  and  Fig.  12  for  13  percent 
slip.  The  computed  and  measured  load,  torque,  and  drawbar  pull 


S! 

5 

.? 


i 


i. 


25 


INTERFACE  STRESSES  (LBS/SQ  FT) 


Fig.  10  Comparison  of  Measured  and  Computed  Interface 

Stresses  for  a  1. 15  ft  Diameter,  0. 36  ft  Wide  Driven 
Wheel  in  Sand  (33%  Slip) .  Run  #102 . 
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INTERFACE  STRESSES  (LBS/SQ  FT) 


Fig.  11  Comparison  of  Measured  and  Computed  Interface  Stresses 
for  a  1. 15  ft  Diameter,  0,  36  ft  Wide  Driven  Wheel  in 
Sand  (22%  Slip)  •  Run  #105. 


Interface  Stress  (Ib/sq  ft) 


values  are  shown  in  Table  1.  The  constants  In  the  slip  relation, 
Eq.  (6),  were  determined  from  the  actual  wheel  performance  tests 

as  1  »  0.098  and  K  ■  0.76.  Further  research  is  needed  to  de- 

o 

t ermine  whether  these  constants  are  generally  valid  for  cohesion- 
less  soils  and  other  wheel  dimensions. 

4.2  Dense  sand.  In  this  test  series  the  sand  was  compacted 
by  vibration.  For  this  condition  the  friction  angle  was  found  to 
be  4>  «  44°  .  Cone  penetrometer  tests  indicated  that  wheel  passage 
actually  loosened  the  sand.  Since  this  loosening  is  the  end  re¬ 
sult  of  the  shearing  process,  the  same  friction  angle  was  used  for 
the  rear  field  as  in  the  front.  Interface  stresses  obtained  in  a 
typical  test  are  shown  in  Fig.  13  and  pertinent  data  are  shewn  in 
Table  1.  The  slip  for  the  dense  sand  was  computed  by  the  con¬ 
stants  determined  for  loose  sand;  the  constants  appear  to  be  valid 
for  both  cases. 

4.3  Loam  (light  clay).  For  this  test  series  the  loam  bed 
was  prepared  at  about  16  percent  moisture  content  and  lightly 
rolled  to  provide  a  uniform  surface.  Strength  properties  of  this 
soil  were  also  determined  by  triaxial  tests  performed  at  a  loading 
rate  comparable  to  that  in  the  wheel  performance  tests.  Cone 
penetrometer  tests  were  conducted  before  and  after  the  wheel  pas¬ 
sage  and  also  in  small  soil  bins  where  placement  and  moisture  was 
well  controlled.  These  tests  served  the  purpose  of  correlating 
triaxial  test  results  and  actual  conditions  in  the  test  bed. 

Figure  14  shows  actual  measurements  from  three  runs  performed 

* 

under  the  same  conditions  and  the  computed  stresses  for  C  «* 

130  lbs/sq  ft,  and  4>  »  23°  in  the  front  field,  and  C  • 

195  lbs/sq  ft,  and  <t>  -  20.7°  in  the  rear  field.  The  interface 
stress  measurements  for  the  three  runs  are  almost  identical  in¬ 
dicating  the  reproducibility  of  the  results.  The  computed  stresses 
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COMPARISON  OF  MEASURED  AND  COUNTED  WHEEL  PERFORMANCE  PARAMETERS 


INTERFACE  STRESSES  (LB/SQ  FT) 
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Fig,  15  Comparison  of  Measured  and  Computed  Interface 

Stresses  for  an8*in  Diameter,  2 -in  Wide  Driven  Wheel 
in  Jones  Beach  Sand-Bin  Sloped  4°  (17%  Slip).  Run  #110602, 
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Fig.  16  ComparrLon  of  Measured  and  Computed  Interface  Stresses 
for  an8«Jtn  Diameter,  2 -in  Wide  Driven  Wheel  to  Jones 
Beach  SaM  -  Bin  Level  (17%  Slip).  Run  No.  110102 


duplicate  the  measured  ones  quite  closely,  although  in  the  ex* 
periments  the  stresses  never  rise  instantaneously,  as  predicted 
by  theory.  The  slip  coefficients  established  for  the  loam  are 
j  •  0.089  and  K  *  0.09,  substantially  different  from  those 
for  sand.  Computed  and  measured  wheel  performance  parameters  are 
shown  in  Table  1. 

4.4  Another  series  of  validation  tests  using  a  smaller 
sized  bin  and  wheel  was  conducted  at  Grumman.  In  this  test 
series,  an  8-inch  diameter,  2-inch  wide,  wheel  was  instru¬ 
mented  to  measure  interface  normal  and  shear  stresses  at  the 
centerline  and  then  at  the  edge  of  the  wheel.  The  bed  material 
was  medium  dense  Jones  Beach  sand  whose  properties  had  been  pre¬ 
viously  determined  and  reported  (Refs.  7  and  9).  The  tests  were 
performed  for  various  loads  at  slips  ranging  from  9  percent 
drag  to  39  percent  slip.  Friction  angle  in  the  forward  zone 
was  36  degrees  and  in  the  rear  zone  41  degrees  based  on  re¬ 
sults  of  triaxial  and  cone  penetrometer  tests.  The  bin  at  Grum¬ 
man  possesses  the  unique  capability  of  providing  a  slope  up  which 
the  wheel  may  be  driven  or  pulled.  Therefore,  the  effects  of 
slope  on  wheel  performance  characteristics  can  be  studied.  Fig¬ 
ures  15  and  16  show  the  measured  and  computed  distribution  of 
average  normal  and  shear  stresses  for  the  level  and  sloped 
(4  degree)  cases,  respectively.  Measured  slip  in  both  cases 
was  17  percent.  The  computed  and  measured  load,  torque,  and 
drawbar  pull  values  are  shown  in  Table  1. 

5.  EXTENSION  OF  THEORY  FOR  THE  PREDICTION  OF  TIRE  PERFORMANCE 

Rigid  wheel  performance  closely  approximates  that  of  tires  if 
the  stiffness  of  the  tire  is  great  relative  to  the  soil.  For  less 


stiff  tires  it  is  necessary  to  take  the  effect  of  tire  deformation 
on  soil  reaction  into  account. 

Tire  deformation  affects  soil- tire  interaction  in  two  ways: 
it  changes  the  geometry  of  the  soil-tire  interface,  and  it  re¬ 
lieves  the  stresses  that  would  develop  in  the  soil  if  the  inter¬ 
face  were  undef ormable . 

A  soil-tire  model  that  allows  the  consideration  of  all  essen¬ 
tial  factors  affecting  tire  performance  is  shown  in  Fig.  17.  The 
tire  is  assumed  to  have  a  constant  width  both  in  the  undeformed 
and  in  the  deformed  state.  The  stresses  across  the  tire  width  are 
assumed  to  be  uniform  so  that  the  soil-tire  interaction  problem 
can  be  treated  as  two  dimensional.  Tire  deformation  is  repre¬ 
sented  by  the  shape  of  the  tire  in  the  center  plane  in  the  direc¬ 
tion  of  travel.  Tire  shape  is  assumed  to  be  the  same  in  all  par¬ 
allel  planes.  The  deformation  is  assumed  to  consist  of  two  curvi¬ 
linear  segments  separated  by  a  linear  or  flat  section.  It  is  as¬ 
sumed  that  the  tire  starts  to  deform  an  angle  a'  ahead  of  the 
entry  angle  (a  ) ,  and  reaches  its  original  form  an  angle  a" 
past  the  exit  angle  a^.  In  the  front  curvilinear  segment,  the 
radii  decrease  according  to  a  logarithmic  relationship. 

The  interface  stresses  between  the  entry  arr-  ?,  ag  and  a^, 
as  well  as  those  between  the  exit  angle,  ar  and  a^,  are  as¬ 
sumed  to  be  controlled  by  failure  conditions  in  the  soil.  Out¬ 
lines  of  the  respective  slip  line  fields  are  shown  by  dashed  lines 
in  Fig.  17.  The  geometry  of  these  failure  zones  and  the  associ¬ 
ated  stresses  are  computed  the  same  way  as  for  rigid  wheels.  In 
the  flat  portion  of  the  tire  the  soil  pressure  equals  P^,  a 
limit  pressure  characteristic  of  the  tire  and  inflation  pressure. 
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DEFLECTION  OF  TIRE  STARTS 


LU 

z 

o 


Fig.  17  Soil-Tire  Model 


Preliminary  work  with  this  model  indicates  that  experimental 
results  can  be  duplicated  with  the  model  reasonably  well  (Ref.  10) 
We  plan  to  incorporate  this  model  in  our  computer  program  so  that 
our  method  of  tractive  force  prediction  would  be  applicable  to 
both  rigid  wheels  and  tires. 


6.  SUMMARY 

The  experimental  ..esults  obtained  for  a  variety  of  soil  con¬ 
ditions  indicate  that  the  concept  of  soil-wheel  interaction  pre¬ 
sented  in  Section  3  is  valid,  ar.d  interface  stresses  developing  be¬ 
neath  rigid  wheels  are  controlled  by  soil  failure  conditions.  The 
numerical  integration  of  the  differential  equations  of  plasticity 
for  soils  is  the  basic  tool  used  in  determining  interface  stresses 
for  the  boundary  conditions  imposed  by  the  wheel.  Application  of 
these  principles  to  the  determination  of  wheel  performance  parame¬ 
ters  results  in  a  predictive  method  directly  related  to  soil 
strength  parameters. 

In  the  development  of  the  theory  of  soil-wheel  interaction  it 
was  necessary  to  make  certain  assumptions  as  to  the  development  of 
the  interface  friction  and  slip.  The  interface  friction  angle, 

6,  was  assumed  to  be  uniform  over  the  perimeter  of  the  wheel. 
Although  in  most  of  the  experiments  this  assumption  was  found  to 
be  approximately  correct,  in  some  experiments  on  clay  the  inter¬ 
face  friction  angle  was  found  to  increase  monotonically  from  the 
entry  angle  toward  the  rear  angle.  Although  such  a  behavior  was 
found  not  to  affect  the  wheel  performance  parameter  prediction 
excessively,  further  research  could  clarify  the  cause  of  such  be¬ 
havior  and  allow  appropriate  modification  of  this  assumption. 
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The  other  area  where  further  experimental  information  would 
be  most  useful  is  the  shear  stress-slip  relationship.  It  would 
be  desirable  to  establish  appropriate  constants  for  the  various 
types  of  soils  as  a  function  of  strength  properties  so  that  such 
relationships  could  be  used  in  predictive  methods. 
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